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Abstract 

An effect of overlapping of exhausted regions around droplets of 
a new phase is analyzed. Several approximation to take overlapping 
into account are suggested. The small parameter responsible for con- 
vergence is extracted. The multi particle effects of overlapping are 
described. 

Introduction 

Kinetics of the first order phase transitions includes many aspects of evolu- 
tion. One of such aspects is the exhaustion of the volume where the embryos 
can appear. Since A.N.Kolmogorov pQ has solved a problem of exhaustion 
for the case of crystallization this problem is considered to be analyzed com- 
pletely. 

Nevertheless in the case of nucleation of a supersaturated vapor into a 
liquid state of droplets kinetics of nucleation is governed by another mecha- 
nisms which depends on the Knudsen number characterizing the nucleating 
system. Under the small Knudsen numbers the kinetics of nucleation was 
described in [2]. Under the big Knudsen numbers kinetics resembles kinet- 
ics of crystallization but with some modifications. These modifications were 
described in [3]. Here the precise solution of this problem is absent. The 
source of difficulties lies in a very specific behavior of exhausted regions. In 
this paper the further investigation of this problem is presented. 
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1 Direct regions 



The rate of nucleation J as it is shown in |3J is proportional to the volume 
Vf ree of the non-exhausted region 

J ~ Vf ree 

The total volume of the system V to t is supposed to be the unit one. Our task 
is to determine Vf ree as accurate as possible. 
Certainly, 

Vf r ee 1 V e xh 

where V ex h is the volume of the region exhausted by already existed droplets. 

Around the formed droplet there will be an exhausted region formed by 
the vapor consumption by this very droplet. This region will be called as a 
"direct region". 

For V ex h one can write the following representation: 

v exh = f J(t r )v4t,t')dt' 

Jo 

Here t = is the initial moment of time, t is the current moment of time. 
The value V c i(t,t f ) is the (mean) volume which is exhausted at the moment 
t from Vf r ee by the embryos formed at t'. For the last value one can suggest 
several approximations. 

1.1 Kolmogorov's result 

In crystallization Kolmogorov justified the following result 

V cl (t,t') = V (t,t') 

Here Vo (£,£') is the volume of the region which is exhausted at t elsewhere 
by embryos formed at t' . 
In crystallization 

v^it-t'Y 

where d is the dimensionality of space. The linear size R ~ V} /d of the direct 
region grows linearly in time. This is the crucial point in justification of the 
Kolmogorov's result. 



2 



In general situation of nucleation asymptotically at t — t' — > oo one can 
see that 

0(t-t') d/2 < V <O(t-t') d 
It means that the linear size grows as 

0(t-t') 1 / 2 <R<0(t-t r ) 

So, the non-linearity is not too strong and this will be the source of several 
approximations considered below. 

In any situation the dependence R(t — t') is known. Ordinary R ~ (t — 
t'f/ 2 . 



1.2 Stochastic approximations 

In stochastic approximation the following approximation 

V d (t,t') = V free (t)V (t,t') 

is taken. 

This approximation implies stochastic overlapping of volumes Vq coming 
from different droplets. Why for the volume V depending oh t and t' the 
overlapping is taken at t? More reasonable is to take the half-sum 

V d (t,t>) = {Vf - {t) + 2 Vfreeif) ] Vo(t,t') 

or 

V d (t,f) = jW free (t")dt"V (t,f) 

The approximate linearity of R(t — t') leads to the approximate similarity 
of all approaches. 



2 Approximation of pair interaction 

Consider the auxiliary problem which is formulated as following: 

Consider the droplet born at t" . The question is what will be the part of 

Vq occupied by direct regions from other droplets. 

One can approximately suppose that outside direct region of this droplet 

the rate of nucleation is not perturbed and inside the direct region the rate 

of nucleation is zero. 
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Then the problem is purely geometrical one. One has simply to calculate 
the overlapping of spheres. The value V over (r",t",r',t') is the volume of 
overlapping for the droplet born at t" in the point r" and the droplet born 
at t' in the point r'. It is calculated purely geometrically. 

Also purely geometrically one can calculate the boundary r' b of the direct 
region and the time £& when the given r" is attained by the boundary of this 
direct region. 

As the result one can easily calculate the total volume of overlapping by 
integration 



Here t b is a minimum of t and t b , Sd is the surface square of a sphere of radius 
r in a space with dimension d. 

The problem is that the supposition that outside direct region of the 
droplet the rate of nucleation is not perturbed and inside the direct region 
the rate of nucleation is zero is no more than an approximation. 

Another approximation is that the rate of nucleation outside is J(t') taken 
from another approximation, for example, from the previous one. Then 



Unfortunately, the rate of nucleation J{t') strongly depends here on r. 
To see this we can analyze the region (it will be called "an inverse region") 
where the region is forbidden to have the free point r, t' . This region can 
cross the direct region for the droplet born at 0, t". 

It means that in the outside region near r' b the rate of nucleation is higher 
than far from r b 

This correlation is missed here. 

3 Approximations for estimates 

The effective approach is to note that d ^> 1 and use this fact. For example 
for d = 3 the droplets appeared until the first quarter of the nucleation period 
in decay conditions occupy the small part of the whole volume and later all 
evolution is governed namely by these droplets. This fact allows to construct 
approximations as it is shown in [I]. 





4 



4 Inverse region 



Inverse region for the droplet born at r',t' is the geometrical place of points 
where the formation is forbidden to have the possibility to form droplet at 
r',t'. Certainly, this region is situated symmetrically to the direct region of 
the same droplet. 

In the case of the Kolmogorov's solution R = A(t — t') with some constant 
A. Then any inverse region can not cross any direct region. 

Having required that at the bottom of inverse region there is no formation 
of droplets one can see that there is no traces of direct regions in the whole 
inverse region. Then the Kolmogorov's approximation really takes place. 

Consider now the general situation. Let P-(t) be a probability that in 
the point r = at t there is a free point. For P-(t) one can write 

P_(t) = f dt' P dr'S d PJt',r') 
Jo Jo 

where P-(t',r') is the probability that the point r',t' is free. 

Let r',t' be inside th inverse region formed at r = 0, t. For P-(t',r') one 
has to construct the inverse region. This region can go outside the inverse 
region of r = 0, t. This is the difficulty For P_(t',r') one can then write 

P_(t', r') = f dt"[ f d?'P_(t", r") + f dP>P_(t", r")} 

" ^ ^exhausted ^ ^non— exhausted 

Here V exhausted is the inverse region based on r', t' which lies inside the inverse 
region based on 0, t and V non - exhausted is the inverse region based on r', t' which 
lies outside the inverse region based on 0, t. 

The last separation is artificial but it will be essential in further approx- 
imations. 

The integral equation is the linear one. Moreover, in dimension d = 3 
there exists a cylinder symmetry. 

The suitable way is an iteration solution based on the smallness of over- 
lapping 

P-(i)(t', r') = f dt"[ f dP'P^it", r")+ f dJP-wW, r")] 

"^0 ^Vexhausted ^ ^non — exhausted 

The initial approximation is evident: P = 1 in V ex hausted, P = Pmean in 
Vnon-exhausted- Here P m ean is some mean rate of nucleation calculated by 
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other approximations. The smallness of overlapping ensures the convergence 
of iterations. 

This consideration is also some approximation. Here we do not take into 
account the possibility to born the third droplet with its own inverse region 
which causes the change of P_. 

5 Approximate simplification 

If one approximates J by initial intensity Jo in V ex hausted and by the mean 
intensity J in V non _ ex h auste d then one can see a very simple expression 

, T ) / dt [JoVexhausted ~l~ -J^non— exhausted] 

Jo 

The last approximation is both simple and accurate one. 

Earlier the effects of interaction (overlapping) of three and more regions 
were not taken into account. Since the effects of overlapping are small, the 
influence of effective interaction of several particles are small. 

6 One limit case 

In order to construct the procedure taking into account effects of overlapping 
(interference) of exhausted regions initiated by different droplets one has to 
outline the parameter responsible for the relative smallness of the many- 
droplets effects. Only then the effects under consideration will have the 
order of corrections and can be considered as some perturbations. 

The small parameter e of the theory will be, thus, a ratio of the square 
of exhausted region with overlapping to the the square of exhausted region 
without overlapping. We expect that e is seriously less than 1. 

Ordinary the growth of the linear size of the exhausted region can be well 
approximated by a power-like dependence 

R(t,t') =A(t-t') a 

with parameters A and a. The case a — >• oo corresponds to the maximum of 
e. 

Fortunately, this case allows analytical solution. 
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The case a = oo corresponds to instantaneous creation of the exhausted 
region of some radius Rq. Then the probability P_(r, t) to be the point r free 
at the moment t equals to 



P_ (r, t) = I dr' f dt'P- (/, t')I 
Js Rn Jo 



where the space integral is taken over the sphere of the radius P with the 
center r and I is the kinetic factor, i.e. the free intensity of the droplets 
formation. 

The crucial point of consideration is that here P__ can not depend on r' . 
Then ^ 

P-(r,t) = 4ttR 3 [ dt'P_(r,t')I /3 
Jo 

The differentiation gives 

^l^ = 4 7 rP 3 P_(r ) t)/ /3 

The last differential equation has an evident solution 

P_(r,t) = exp(-4nRlt/3) 

where initial moment of time is taken as t = 0. 

Then it is easy to calculate the total number of exhausted regions 

M tot = dtP-(r,t)=J dtexp(-A7rR 3 t/3) = 3/(4ttP 3 ) 

The absence of overlapping corresponds to the linearized version, i.e. 

r 3/(47rR3) 



Mi 



Itn 



== / dt(l - 4ttP^/ 3 ) = 3/(8vrP||) 

Jo 



So, M to t = 2Mn n . Since the squares of all exhausted regions are 47tPq/3, the 
last ratio will be the power of overlapping, i.e. 

e = (Mtot - M lin )/M lin 

So, in any situation e and can be considered as a small parameter. 
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7 The scheme of calculation 



Now it is worth to present a concrete effective scheme of construction of the 
multi-particle effects of overlapping. 

The first step is the formal one - to form the first droplet at r',t'. 

The rate of nucleation inside the exhausted region of the first droplet 
(ERi) is zero. The question is to find the rate of nucleation outside ER\. 

Then it is possible to choose two alternative possibilities to go forward. 
The first one is to choose the rate of nucleation outside ER\ to be the ideal 
value I . This approach is very close to the direct simulation of the process. 
Unfortunately, in this approach to get a suitable effect it is necessary to take 
into account many neighbor droplets. For example, in the plane case the 
number of neighbors has to be greater than six. In the free-dimensional case 
this number is greater than twelve. 

The second possibility is more effective. Instead of I outside the ER 1 
one can use the average rate of nucleation /. One can define the average rate 
of nucleation I v over the volume V as 

1 



= y J drl(r,t) 



The value / is the limit value 



' v 



Namely this value is chosen as the rate of nucleation outside ER X . 

Corrections to the last value appear because near the boundary of ER\ 
some droplets which might affect on a point outside ER\ had been born in 
ER\ where formation is forbidden. 

The rate of this value has an order of relative volume of intersection of 
ER\ and inverse region IR of the second droplet, i.e. IR2. More precise the 
probability P2 to form the droplet at r", t" when the first droplet is formed 
at r' t 1 

Here Vol means the volume. The value Vol(ER 1 nlR 2 ) is the volume of the 
intersection of ERi and IR 2 , etc. 

The volumes standing here are the simple known algebraic functions, the 
calculation is evident. 
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The effectiveness of this procedure is based on the relative smallness of 
the value 

Vol(ER 1 n IR 2 ) 



a = 



Vol(IR 2 ) 

The average probabilities are given by 

Since the probability Pi(r',t') does not depend on r' the integral over r' 
is not more than a formal object. 

It is more convenient to consider the relative deviation of P 2 from the 



base value 7- 



AP 2 (r»X;r>t>) = V ° l ( ER i nIR ^ 



AP 2 (r",t";t') = v J v dr 



Vol{IR 2 ) 
1 r ,Vol(ER 1 fMR 2 ) 



Vol{IR 2 ) 



v ' ; Jo V Jv Vol(IR 2 ) 
The analogous constructions can be made for the three ER overlapping. 
The relative deviation in probability 5 P 3 {r"' , t" 1 '; r" , t" '; r' , t') connected 
with simultaneous overlapping of three volumes is given by 

5P 3 {r ,t ;r,t;r,t) = Vol(IR 3 ) P2{r,t;r,t) 

Here there is the deviation in probability to find the third particle born at 
r'", t'" when the first particle was born at r', t' and the second particle was 
born at r", t" . 

The relative deviation in probability 5P 3 (r"', t'"; r", t"; r', t') simply to find 
three particles born in positions r'", t" r ; r", t"; r', t' is given by 

6P s (r'", t>"; r", t"; r>, f) = ™ (( ^ ( fff ° '^ PtW, t"; r> , f) 

The knowledge of P\,P 2 , PS solves the problem to find the rate of nucle- 
ation is one-particle, two-particle and three-particle approximation. 

Certainly ER(t', r') as a function of t' — t instead of t coincides with 
IR(t',r') 
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8 Volumes of overlapping 

In two particle approximation the volume of overlapping V 2 can be calculated 
in a very simple way 

V 2 (r'.t') = J dt" J dr"Vol{ER l nER 2 )P 2 {r",t";r',t') 

Here the argument r' is not important because the system is uniform. 

The partial effective volume in the two-particles approximation looks like 

V eff = Vol(ER 1 )-±V 2 (r',t') 

In the three-particles approximation the volume of overlapping of three 
volumes is given by 

V 123 = J dt" J dr" J t!" J dr'"Vol(ER 1 nER 2 nER 3 P 2 (r",t";r'.t')P 3 (r'".t'";r".t 

It is necessary also to introduce the volumes of overlapping of ER of two 
particles in the presence of the third one. The indication of a particle which 
does not have ER in the two-particle overlapping stands in brackets. Then 

Vl2{3) = J dt" j dr" J dt'" J dr"'Vol(ER 1 nER 2 )P 2 (r".t";r',t')P 3 (r"',t"';r",t"; 

The volume of two-particles overlapping without the volume of the three- 
particles overlapping is approximately given by 

\> 12(3) = J dt" j dr" J dt'" j dr"'Vol((ER l n ER 2 ) - (ER 1 n ER 2 n ER 3 )) 

P 2 (r".t";r',t')P3(r"',t"';r",t";r',t') 
Analogous value is given by 

V 1{2)3 = j dt" j dr" j dt'" J dr"'Vol((ER 1 n £i? 3 ) - (ER l n ER 2 n ER 3 )) 

P 2 (r".t";r',t')P 3 (r"',t"'-r",t";r',t') 
The effective volume in three particles approximation is given by 

V eff = VoliER,) - V 12{3) /2 - V 1{2)3 /2 - 2V 123 /3 

One has to stress that all volumes indicated as Vol can be elementary 
calculated. 
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